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We study the complementarity between the indirect detection of dark matter with 7-rays in 
H.E.S.S. and the supersymmetry searches with ATLAS at the Large Hadron Collider in the Fo- 
cus Point region within the mSUGRA framework. The sensitivity of the central telescope of the 
H.E.S.S. II experiment with an energy threshold of ~ 20 GeV is investigated. We show that the de- 
tection of 7-ray fluxes of O(10 -12 ) cm _2 s _1 with H.E.S.S. II covers a substantial part of the Focus 
Point region which may be more difficult for LHC experiments. Despite the presence of multi- 
TeV scalars, we show that LHC will be sensitive to a complementary part of this region through 
three body NLSP leptonic decays. This interesting complementarity between H.E.S.S. II and LHC 
searches is further highlighted in terms of the gluino mass and the two lightest neutralino mass 
difference. 
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I. INTRODUCTION 



A preferred particle physics candidate for the Dark 
Matter (DM) component in the Universe is an electri- 
cally neutral Lightest Supersymmetric Particle (LSP), 
which in various supersymmetry (SUSY) breaking sce- 
narii is the lightest neutralino x GL 01- The se ^~ 
annihilations of neutralinos in halos of galaxies would 
produce Standard Model particles: 7, v, charged lep- 
tons and hadrons. The ongoing and forthcoming ground- 
based (H.E.S.S. 3, MAGIC [|, VERITAS [| and space 
(AMS @, GLAST Q, PAMELA 0]) experiments de- 
tecting galactic cosmic rays will open large new windows 
for the DM annihilation signal detection. In particular, 
phase II of the H.E.S.S. experiment will provide an im- 
proved sensitivity due to a combination of a large effec- 
tive area, low energy threshold, and an excellent angular 
resolution necessary for the efficient hadron background 
suppression. 

Simultaneously to the new generation of the cosmic ray 
detectors which will provide outstanding data in the 
dark matter detection domain, the Large Hadron Col- 
lider (LHC) start-up in 2008 will explore physics in the 
TeV energy range and thus probe the origin of elec- 
troweak symmetry breaking. It will also introduce a new 
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era of physics and particle searches beyond the Standard 
Model, and in particular ought to discover the new par- 
ticles predicted by supersymmetry or other new states 
such as those foreseen by models with extra spatial di- 
mensions. Under the hypothesis that dark matter is com- 
posed of a SUSY (or Kaluza-Klein) particle, a relation- 
ship between the astrophysical signal observation and the 
properties of the new particles found at LHC can be es- 
tablished. The information from LHC measurements will 
also allow us to discard various hypotheses on the nature 
of the DM [{|. On the other hand, in case of a signal 
detection which is not correlated with an astrophysical 
source, H.E.S.S. would be able to investigate the DM dis- 
tribution profile in halos of galaxies independently of the 
particle physics uncertainties. Accurate reconstruction 
of the morphology of the sources could help to discard 
specific halo models. 

The aim of this paper is to study some experimental 
and phenomenological aspects of the complementarity 
between searches of SUSY particles in the ATLAS ex- 
periment and the detection potential of the H.E.S.S. tele- 
scope in phase II, in the case of the so-called Focus Point 
(FP) region in the mSUGRA framework of the Mini- 
mal Supersymmetric Model (MSSM). The most realis- 
tic Monte Carlo simulations for both experiments will 
provide input to this study. As the measured variables 
in ATLAS are related to various mass differences of the 
SUSY states, the H.E.S.S. II sensitivity will be investi- 
gated with respect to these variables. In particular, we 
will investigate the dependence of the gamma-ray fluxes 
and the annihilation cross-section on the gluino mass and 
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the mass difference of the two lightest neutralinos. De- 
tailed studies in the literature have investigated the com- 
plementarity of accelerator measurements at LHC/ILC 
and astrophysical observations, e.g. [l(| [HI, to unveil 
the nature of Dark Matter. Here, we focus on the com- 
plementary constraints in the FP region with the forth- 
coming atmospheric Cherenkov telescope H.E.S.S. II and 
LHC searches using realistic performances for both ex- 
periments. This FP region is also of great interest for 
searches via neutrino telescopes [12| . 

The rest of the paper is organized as follows: some key 
points of the supersymmetric focus point framework and 
related phenomenological features are given in section 2. 
Section 3 describes various issues related to 7-ray signals 
expected in Cherenkov telescopes and the detection sen- 
sitivity of the upcoming phase II of the H.E.S.S. exper- 
iment. The constraints from the 7-ray flux senstivity to 
the parameter space in the focus point region are studied 
in section 4. MSSM spectrum measurements at the LHC 
for this parameter space region are discussed in section 
5. Section 6 is devoted to the potential complementarity 
between H.E.S.S. II and ATLAS measurements. Conclu- 
sions and perspectives are given in section 7. 



II. SUPERSYMMETRIC FRAMEWORK AND 
THE FOCUS POINT REGION 

We consider the framework where supersymmetry 
breaking in a hidden sector is mediated to the visi- 
ble world through gravitational interactions between the 
very heavy hidden fields and the MSSM sector. Although 
little is known about the actual dynamics responsible for 
this breaking and its mediation mechanisms, it is reason- 
able to assume that the gross features of the low energy 
spectrum can be encapsulated in a universality assump- 
tion (and possible deviations from it) for the effective soft 
supersymmetry breaking parameters at a scale close to 
the Planck or the GUT scale. The more detailed features 
of this spectrum are then determined by quantum effects 
obtained through the running down to the electroweak 
scale. It is important to keep in mind that the theo- 
retical uncertainties at the high (SUSY breaking) scale 
are different in nature from those involving the running 
to the low (electroweak) scale. The former cannot be 
reduced without a better knowledge of the underlying 
new physics, while the latter are in principle reducible 
through the improved calculational theoretical tools such 
as the Renormalization Group Evolution (RGE) codes 
and spectrum calculators, where for instance scale depen- 
dencies, threshold effects, higher order corrections, etc... 
are standardized. In practice, the effects of these two 
types of uncertainties are not easy to disentangle when 
one tries to determine in a model-independent bottom- 
up approach the fundamental parameters of the model, 
ultimately starting from experimental data. Moreover, 
in a collider environment such as the LHC one would 
presumably aim first at typical SUSY discovery events, 



then at the determination of mass differences of vari- 
ous SUSY particles which are kinematically accessible. 
With this respect, resorting to detailed predictions in a 
top-down approach should be taken only as a guide for 
the possible patterns of SUSY mass spectra, production 
rates, branching ratios, etc., as well as the complemen- 
tary requirements for a good dark matter candidate (relic 
density, predictions for direct and indirect dark matter 
searches). 

In this paper we focus on MSSM spectrum patterns 
where all the scalar quarks and leptons are in the multi- 
TeV range and thus out of direct reach at the LHC, while 
the lighter neutralinos and charginos and of course the 
lightest Higgs remain accessible. Such patterns are moti- 
vated by the so-called focus point scenarios [H[ where 
it is noted within mSUGRA that for tan/3 > 5 and 
Wi/2> A) ^ 1 TeV, mo could be very large and still al- 
low electroweak symmetry breaking and the right Z bo- 
son mass to occur without a large fine-tuning between 
the order parameter m 2 H (< 0) and fJ?(> 0) at the elec- 
troweak scale. This feature is due to the quantum ef- 
fects which push the running of m 2 H to an essentially 
unique value, a focus point, at some given low scale, ir- 
respective of its initial value at very high scales. It so 
happens that for the above range of parameters this fo- 
cus point occurs at a scale Qo of order the electroweak 
scale, and where |7Tijy u (Qo)| is also of the same order, 
together with H 2 {Qo)(~ Tn% (Qo)) due to electroweak 
symmetry breaking (EWSB). In this case, the smallness 
of the fi parameter leading to a non negligible higgsino 
component for the lightest neutralino has immediate con- 
sequences on the dark matter issues [H], EH , reducing the 
neutralino relic density and increasing the reaction rates 
and fluxes respectively for direct and indirect dark mat- 
ter detection. These distinctive features, together with 
other ones due to the heaviness of the sfermion sector 
(such as the suppression of flavor changing neutral cur- 
rents and the electric dipole moments of the electron and 
the neutron) makes this scenario with low fine-tuning 
very interesting to study. However, the low fine-tuning 
in the EWSB condition seems to be traded for a very 
large sensitivity to the top yukawa coupling and thus 
to the phy sical top mass, as was shown in various stud- 
ies in [l6| (see also HH). The practical consequence of 
this sensitivity is the important dispersion in the out- 
put of different spectrum calculator and RGE codes and 
thus in the predictions of the various experimental ob- 
servables, not to mention numerical convergence issues 
as pointed out in fl5f through a comparison of the two 
codes ISAJET 7.69~g3 and SUSPECT 2.34 0. These 
features hint at the need for further standardisation of 
the various codes through the inclusion of improved cal- 
culations at similar levels; but we anticipate that this 
would probably not be enough to improve substantially 
the situation, due to the intrinsic high sensitivity on in- 
put parameters in the focus point regions. For that, one 
would require a better theoretical understanding of this 
high sensitivity (meaning of the dependence of the focus 
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point location on the running scale, etc). Since in the 
present study we are mainly interested in the potential 
of the experimental complementarity between H.E.S.S. 
phase II and ATLAS for the challenging multi-TeV spec- 
trum pattern, we will ignore the above theoretical uncer- 
tainties and consider the focus-point-like regions sticking 
for convenience to a given code, namely ISAJET |73] to 
generate the MSSM mass spectrum and couplings. We 
have however performed cross-checks with the SUSPECT 
code [H[ and found reasonable agreement for the overall 
focus point region (see also section 4 for further discus- 
sions) . The mSUGRA parameter space is defined by the 
following set of four parameters and a sign[58(: 



m , m 1/2 , A Q , tan/3, sign{n) 



(1) 



with mo the common scalar mass, mi« the common 
gaugino mass, A$ the trilinear coupling, sign(fi) the sign 
of the higgsino mass and tan /3 the Higgs vacuum expec- 
tation values (VEV) ratio. In the sequel of the paper, 
we will perform detailed scans of the parameter space in 
the focus point region using both ISAJET 7.69 interfaced 
with DarkSUSY 4.1 Fiji, and ISAJET 7.71. 



III. INDIRECT DETECTION OF DARK 
MATTER 

A. Galactic sources of dark matter halos 

The choice of the astrophysical targets to be ob- 
served is of main importance for the dark matter sig- 
nal searches. A mong the most commonly s ugg ested 
sources [2(], HH, [H, [23j |. various studies such as [2Cj have 
shown that the Galactic Center could be an interesting 
object under the assumption that the DM density halo 
presents a strong density increase towards its inner re- 
gion. The dark matter density profile extrapolated to 
the central region is currently parameterized following 
the cored spherical models as proposed by [24[ or as halos 
with a cusp as proposed by [25| and [26[ . The cuspiness 
of the halo is not yet confirmed by the rotation curves of 
the stars [27j and is in contradiction with a high value of 
the microlensing optical depth in the central region [28[ 
which cannot be due to the diffuse DM but rather to 
the presence of faint stars and brown dwarfs. The other 
argument against the presence of a large density of the 
DM in the Galactic Center is its impact on the large 
galactic bar rotation which could be modified by the dy- 
namical friction with dark matter medium [29j. However, 
microlensing optical depth and dynamical friction argu- 
ments are slightly controversial [30j and a cusp at the 
Galactic Center is not yet ruled out. It has to be under- 
lined that the Galactic Center region is rich in potential 
7-ray emitters such as supernova remnants, the super- 
massive black hole Sgr A*, the recently discovered HESS 
J 1745-290 source [31(, and the standard interactions of 
charged cosmic rays producing diffuse 7-ray emission [32j . 



The most promising astrophysical objects with large 
mass-to-luminosity ratio are the Spheroidal Dwarf Galax- 
ies (dSph) of the Local Group, periodically crossing the 
Galactic Plane. The not too distant ones are Draco, 
Sagittarius and Canis Major located respectively at ~ 
80, 24 and 8 kpc from the heliocentric position. The 
presently measured star rotation curves provide the 
mass-to-luminosity ratio in the range of 20 to 100 and 
no final answer has been given about cuspiness or not of 
the dark matter density profile, so the cored density pro- 
file should also be considered when predictions for 7-ray 
fluxes from the neutralino annihilations are performed. 
As the line-of-sight integration of dark matter density is 
mainly driven by the distance value of the object, the Ca- 
nis Major galaxy is certainly the most promising target. 
On the other hand, this galaxy seems to be the most dis- 
rupted one by the tidal forces as it orbits our Milky Way. 
The expected flux ratios between dSphs and the Galac- 
tic Center taken as a reference source are respectively ~ 
0.03, 0.15 and 0.5 for Draco, Sagittarius and Canis Major, 
as calculated for a Navarro-Frenk- White (NFW) cusped 
profile by (HI ■ Despite smaller 7-ray fluxes expected for 
dwarf galaxies, they consist of less complex environments 
compared to the Galactic Center region as pointed above 
and present lower background coming from astrophysical 
sources. The extended TeV emission along the galactic 
plane leads to a diffuse astrophysical background which 
is very challenging to overcome. Dwarf galaxies are cur- 
rently in the scope of the H.E.S.S. experiment and first 
results on Sagittarius have been reported in [34j |. 

In this paper we consider the Galactic Center as a 
benchmark source in order to compare our predictions 
with those provided by other authors. The 7-ray fluxes 
from annihilation of neutralinos in a spherical dark halo 
are obtained from 



d$ 7 (Afi,£: 7 ) 1 (cru) dN^ 



dE~, 



4tt 2 m 2 x dE 1 



xJ(Aft)An (2) 



with 



J W = An 



dn 



An 



p 2 \r 



(s)]ds 



(3) 



Eq. ((5J) is expressed as a product of a particle physics 
term and an astrophysics term. The former contains 
(<rv), the velocity- weighted annihilation cross-section, 
the neutralino mass m x and the gamma-ray annihila- 
tion spectrum dN 1 /dE 1 . The astrophysics part depends 
on p, the radial density profile expressed in terms of 
Sq — 2ssq cos 9 where s is the heliocentric 
distance of a given point along the line of sight (l.o.s) in 
the Galactic halo, So is the heliocentric distance of the 
Galactic Center, and 9 denotes the angle between the di- 
rection of the Galactic Center and the l.o.s. The integral 



over s ranges from to sq cos 9 



R l 



sin 2 6 where 



i?o is the radial extension of the spherical halo. The in- 
tegration is performed along the l.o.s. to the target and 
averaged over the solid angle AVl usually matching the 
angular resolution of the detector. 
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B. The H.E.S.S. experiment 

The High Energy Stereoscopic System (H.E.S.S.) con- 
sists of four Imaging Atmospheric Cherenkov Telescopes 
(IACTs) 0. It is designed to detect very high energy 
(VHE) 7-rays in the energy range from 100 GeV up to 
100 TeV. H.E.S.S. detects Cherenkov light emitted from 
electromagnetic cascades of secondary particles resulting 
from the 7-ray/hadron primary interaction in the up- 
per atmosphere. Cherenkov images of air showers are 
used to deduce the energy and direction of the primary 
particle. At the zenith, the energy threshold of the sys- 
tem is 100 GeV and for point sources, an energy resolu- 
tion of 15% is achieved with stereoscopic measurements 
which determine the height of the maximum shower de- 
velopment. The angular resolution for individual 7-ray 
is better than 0.1° and the point source sensivity for a 
5<T detection reaches 1% of the flux of the Crab nebula 
in 25 hours [35|. The combination of these unequalled 
characteristics allows for detailed studies of high energy 
7 source morphology. 

The future of H.E.S.S. is in an u pgr ade to the existing 
telescope array. H.E.S.S. phase II [36( consists of a very 
large single telescope located at the center of the H.E.S.S. 

I array, expected to start taking data as soon as end of 
2008. The total mirror collection area is about 600 m 2 
and the camera with a 3° field of view is made up of 
2048 photomultiplier tubes. The pixel size of 0.07° will 
provide a better resolved shower image as compared to 
H.E.S.S. I. In the stand-alone operating mode, H.E.S.S. 

II will reach an energy threshold as low as 20 GeV. 



C. Detection sensitivity of ground-based 
Cherenkov telescopes 

The phase II of the ground-based Imaging Atmospheric 
Cherenkov Telescopes with its improved sensitivity at en- 
ergies below 100 GeV and excellent angular resolution, 
may allow exploration of significant portions of the SUSY 
parameter space. Its sensitivity will depend strongly on 
the hadron background suppression procedures that will 
be developed at a few tens of GeV, unlike the space tele- 
scopes where the background is dominated by the diffuse 
photon emission. At present energies of H.E.S.S. I above 
100 GeV threshold, the photon/hadron shower identifi- 
cation allows us to discover weak extended sources in a 
reasonable time of observation below 50 hours. One has 
to keep in mind that the phase II of H.E.S.S. will provide 
new results in the interesting astrophysical and/or new 
physics domains only if the steep hadron (proton, nuclei, 
etc) and cosmic electron spectra can be discarded. Forth- 
coming ACTs (such as the MAGIC II experiment [4j with 
two 17 m diameter telescopes) will investigate energies 
below 100 GeV. However, such ACTs are not expected 
to lower their energy thresholds as much as H.E.S.S. II 
will do, since these thresholds scale roughly linearly with 
the mirrors' size. 



The following observation conditions of a given source 
will determine the effective sensitivity to the SUSY sig- 
nal: 

• the elevation angle as seen by H.E.S.S. on the hori- 
zon will have strong impact on the energy thresh- 
old, the hadron background suppression and energy 
resolution; 

• the number of active telescopes, the zenith angle 
corresponding to a given source observation and 
its angular offset with respect to the center of the 
camera, will determine the acceptance values as a 
function of energy; 

• the source spatial profile and weakness of the signal 
will condition the accepted number of 7 and the 
power of the background suppression. 

The optimal source observation conditions are currently 
fulfilled in H.E.S.S. experiment by the Galactic Center 
and Dwarf Spheroidal Galaxy observations. 

Typical uncertainties on the integrated flux measured 
with H.E.S.S. are of the order of 20% which account 
mainly for systematics due to different signal extraction 
methods. The main contribution to these systematic ef- 
fects come from the uncertainties on the hadron back- 
ground suppression factors, strongly varying with the 
studied energy domain. As far as the low energy part 
of the spectrum is concerned, the data to come from 
H.E.S.S. phase II will allow us to explore the 20 GeV 
energy range for which there is no prior experience. The 
values quoted here are provided by Monte Carlo studies 
relying on the extrapolation from higher energies. The 
impact of the energy resolution on the energy threshold 
determination may be considered as negligible. The fore- 
seen experimental uncertainty is still much less than the 
astrophysical uncertainties (see also section 4). 



D. Sensitivity calculation 

In the case of Cherenkov telescopes where the signal is 
extracted by ON-OFF subtraction, the significance S is 
given by [37| 



S 



V, 



V2Afc 



(4) 



ck 



where iV 7 is the number of detected photons from DM 
annihilations, and the number of background pho- 
tons. 

For each annihilation of neutralino pairs, the differen- 
tial continuum photon spectrum dN^ /dE^ expected from 
SUSY signals can be approximated [38|, S0| : 
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This leading-log approximation is obtained from a quark 
fragmentation model for the total hadron spectrum using 
the Hill spectrum. For more details, see [HI]. On the 
other hand, the number of continuum photons above the 
energy threshold E t h collected by the telescope array can 
be expressed in terms of: 



JV 7 (£ 7 > Eth) = T c 



A eff{ E ~t)-^T dE l 



(6) 



where A e ff(E~ f ) corresponds to the effective area of the 
H.E.S.S. instrument for given energy, zenith angle, etc., 
T b s the observation time, and d$ 7 /d£' 7 the differential 
7-ray flux from DM particle annihilations calculated with 
the formula given in Eq. Using Eq. ((SJ) and Eq. © 
and inserting the expression of d<I> 7 /<iE 7 , after integra- 
tion one can compute the minimum detectable annihila- 
tion rate {av) mm : 



An 



J (AST) ASl 



dN 

A eff( E ~f)-^rdE~, 



(7) 



For the signal under consideration, the major sources of 
background are protons, nuclei and electrons. Since the 
nuclei background is subdominant compared to others, 
this background will be neglected in what follows. 

In the case of the hadronic background, the computa- 
tion of the number of background events to be detected 
is derived using the following expression [2(| : 



had 



dSl 
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(8) 



where thad represents the expected hadronic rejection 
which is in the case of H.E.S.S. II of the order of 80%. 
The contribution from cosmic ray electrons that initiate 
showers, which is not distinguishable from 7-rays, is [2(| 



dSl 
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x 10" 



Eth V 2 - 3 
lGevJ 

[cm~ 2 s _1 sr _1 ] 



TABLE I: The case study parameter region in the Focus point 
regime. The top mass mt is set to 175 GeV (see text for 
discussions on sensitivity to m t .) 



Parameter 



Minimum 



Maximum 



m (GeV) 
mi /a (GeV) 



2200 
150 

sign(n) = + 1 
tan/3 = 10 
A o =0 



5200 
750 



with asymptotically flat velocity dispersion curve, has 
been also considered. Table |TT] reports the value of the 
astrophysical quantity J defined in Eq. ([3]) for the Galac- 
tic Center, in the case of the two aforementioned halo 
profiles for comparison. The resulting range of variation 
of J allows us to estimate the magnitude of the astro- 
physical uncertainties affecting the determination of the 
7-ray flux from neutralino annihilation. Fig. [1] presents 
the 5<T exclusion limit of H.E.S.S. II for the annihila- 
tion cross-section using a given acceptance parametriza- 
tion A e ff(Ej), calculated for a generic Galactic Center 
source with a NFW and a cored DM profile respectively. 
Assuming an observation time of 50 hours and a 20 GeV 
energy threshold, a sensitivity as low as 0(1O -26 ) cm 3 s _1 
for the velocity-weighted annihilation cross-section can 
be achieved in case of a NFW profile. The SUSY mod- 
els from Table U are also shown. In the scanned region, 
neutralino masses range from 50 GeV up to 300 GeV 
and velocity-weighted cross-sections from 10 -29 up to 
~ 10 -25 cm 3 s _1 . The lowest neutralino masses as well as 



TABLE II: Values of J(AQ) for the Galactic Center with a 
NFW and cored profiles respectively [13] for the solid angle 
An = 2 x 10" 5 sr. 





Profile 


J (10 23 GeV 2 cm- 5 ) 




NFW 


270 




Core 


0.7 


(9) 







IV. CONSTRAINTS FROM H.E.S.S. II ON 
GAMMA-RAY FLUXES 

In the framework described in Sec. |H] the SUSY spec- 
trum has been computed for various configurations in 
the Focus Point region. The scanned parameter region 
is defined in Table [U For each set of parameters, the 
integrated 7-ray flux has been derived using Eq. @ in- 
tegrated above the energy threshold. We assume as a 
benchmark halo shape a NFW profile which corresponds 
to a standard distribution for dark matter. A second less 
optimistic halo parametrization, presenting a core profile 



annihilation cross-sections are reached through the mo 
noenergetic lines resulting from loop-induced processes 
such as XX ~> 77, XX lZ and \X -> jh. A large 
fraction of the models obtained in the scan are within 
the reach of H.E.S.S. II in the case of a NFW halo pro- 
file. This opens the possibility to test a fraction of these 
models which satisfy constraints on the CDM relic den- 
sity coming from the Wilkinson Microwave Anisotropy 
Probe (WMAP) [39[. As for comparison, in the case of 
the cored halo profile, the H.E.S.S. II sensitivity reaches 
~ 10 _23 cm 3 s~ 1 and will not be able to constrain this 
parameter space region. 
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The same model set, defined in TableQ] is shown in the 
left-hand side panel of Fig. [3] in the (mo — plane 
for a NFW profile. The higgsino content of the lightest 
neutralino \ implies large annihilation via W + W~ , ZZ 
gauge bosons and substantial amount of 7-rays in the 
final state. For illustration, a typical value of mo as high 
as 3000 GeV is required to obtain a satisfactory value of 
the neutralino relic density for rrit = 175 GeV and TO1/2 
= 300 GeV. Furthermore, as the annihilation channels 
leading to 7-rays control also the thermal relic density, 
we observe a strong (anti)correlation between the latter 
and the 7-ray flux tfhy. As can be seen from left panel 
of Fig. [3l cA 7 spans 5 orders of magnitude, while the relic 
density can vary rapidly over 3 orders of magnitude in the 
considered parameter space region as shown in Fig. [2l 

Only a narrow region can accommodate the WMAP 
constraints on the CDM relic density. In order to ac- 
count for these constraints, we highlight the models 
yielding a neutralino relic density lying in the range 
0.067 < fl x h 2 < 0.156 (corresponding roughly to 5 stan- 
dard deviations). The FP region is extremely sensitive 
to the value of the top mass mt. Indeed, for large mo 
values the steep running of the mn u Higgs doublet mass 
term translates into a high sensitivity of the [i parameter 
to m t through the requirement of radiative electroweak 




m 1/2 (GeV) 

FIG. 2: Neutralino relic density, Q. x h 2 , in the (mo,m 1 /2) plane 
for the FP region defined in Table U with Ao = 0, tan/3 = 10, 
sign(ju) = +1 and m< = 175 GeV. Overlaid is the H.E.S.S. II 
sensitivity for a 5cr detection in 50 hours (black contours). 
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FIG. 1: 5a exclusion limits for H.E.S.S. II on the velocity- 
weighted cross-section (cry) as a function of the neutralino 
mass m x for the Galactic Center with a NFW profile (dashed 
red line). Also indicated is the case of a cored profile (dotted 
green line). The energy threshold is 20 GeV and the obser- 
vation time is 50 hours. mSUGRA models (grey points) from 
the FP scanned region defined in Table U are presented as 
well as those satisfying WMAP constraints on the CDM relic 
density (blue points). 



symmetry breaking. As already mentioned in section 2, 
one can then obtain sustantially different sparticle spec- 
tra when using different RGE codes [4l[ . It is well known 
that the above-mentioned strong sensitivity can lead to 
dramatically different fl x h 2 values for a fixed choice of 
the mSUGRA input parameters mo,mi/2, etc. Thus us- 
ing different codes can possibly require slight readjust- 
ment of the input parameter values to retrieve the ap- 
propriate WMAP allowed region. Nevertheless once a 
scan is performed over a reasonably large domain in the 
(mo,?7ii/2) plane, the effect of using different codes is 
simply to "shift around" a bit the different contours for 
the WMAP allowed ranges, eventually modifying some- 
what their shapes, which has been checked to some ex- 
tent in the present study using the code SuSpect [l|[. 
Thus the discrepancies in the sparticle spectrum in the 
focus region, as obtained from the different publically 
available codes, should not affect appreciably our main 
results, provided of course that the same code is used 
consistently for the whole study. [5^| However, we will oc- 
casionally present the output of two versions of IS A JET, 
7.69 and 7.71, in order to better illustrate the range of 
reliability of the predicted observables. 

The phenomenology of the FP region differs from the 
so-called bulk and co-annihilation regions due to the large 
masses of the sfermions. The scalar sector, including the 
CP even/odd and charged Higgs, lies in the few TeV 
range except for the lightest Higgs boson. In performing 
the scanning of SUGRA parameters, we have checked 
that a number of present experimental and theoretical 
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constraints are fulfilled. Apart from the direct lower lim- 
its on the sparticle masses obtained at LEP or Teva- 
tron [3| which are obviously fulfilled in the range of 
parameters we consider here, there are potentially some 
constraints from virtual supersymmmetric contributions 
to low energy observables. The leading x^t and tH ± 
loop corrections to the branching ratio for radiative b de- 
cays, which can give in principle interesting constraints 
on supersymmetric models, are completely negligible in 
the present case due to the decoupling effects of the very 
heavy sfermion and charged Higgs masses. Similarly, 
supersymmetric contributions [451 ] to = (g^ — 2)/2, 
which would be typically enhanced for large values of 
the combination [i x tan/3 through charginos/sneutrinos 
(and more moderately neutralinos/smuons) one- loop cor- 
rections, are suppressed in our case due to the very large 
sfermion masses and moderate values of tan/3. In par- 
ticular this focus point region would be consistent with 
the comparison between the measurements [46| and the 
standard model theoretical predictions for a M , if the 
hadronic r-decay data (rather than the e + e~ annihila- 
tion into hadrons) are used in the determination of the 
leading hadronic contribution [48| . 

The left panel of Fig. 0] shows the range of gluino 
masses, ro§, as a function of the two lightest neutralino 
mass difference, m„o — m^o, which is a key observable 

' X2 Xi' J 

for ATLAS (see section 5), and the corresponding 7-ray 
flux. The latter observable is indirectly sensitive to the 
neutralino mass difference and somewhat in a model de- 
pendent way. In our mSUGRA constrained scenario, and 
given the LSP mass range and large too values consid- 
ered here, the ZZ and W + W~ annihilation channels 
are open and become actually dominant as compared 
to the fermion pair (//) final state channels. This is 
a combined effect of, on one hand the suppression of 
// final states due to t— channel exchange of very heavy 
sfermions, and on the other hand a substantial higgsino 
component of the LSP enhancing its coupling to Z's and 
W's. A larger m x o — m x o corresponds to heavier neu- 
tralinos (resp. charginos) which are exchanged in the 
t— channel, and thus reduces the ZZ (resp. VF+IU - ) final 
states, implying an overall reduction of the 7-ray fluxes. 
The sensitivity of H.E.S.S. II is displayed as solid black 
contour. Models satisfying WMAP constraints on Q, x h 2 
are also highlighted as dashed purple contours. Gluino 
masses lying in the range 800-1800 GeV are within the 
reach of H.E.S.S. II sensitivity for neutralino mass dif- 
ferences up to 80 GeV. Given the WMAP constraints, 
a significant fraction of cosmologically interesting SUSY 
models will be tested with H.E.S.S. II. 



V. LHC MEASUREMENTS 

The scan described in the previous sections was per- 
formed over a region of mSUGRA space for which rrig < 
1800 GeV. Various studies are available in the litera- 
ture [5(| HH, HH which demonstrate that for an in- 



tegrated luminosity of ~ 100 fb _1 inclusive searches for 
events including multiple jets and/or leptons and missing 
transverse energy would allow the discovery at the LHC 
of gluino production up to masses of ~ 1800 GeV, thus 
covering most of the investigated region. 

Beyond the mere discovery of a signal, the LHC should 
be able to perform measurements of the SUSY spectrum 
which would allow the experiments to put constraints on 
the predicted relic density [l(| H3, [5f| . One point in the 
region addressed here was studied in detail in [4!|, and 
the possible ways of extracting a parameter measurement 
were investigated. It turned out that the cleanest signal 
for parameter measurements would be the study of the 
three-body decay 

X2(3) - (10) 

where the ^2(3) would be produced in the decay g — > 
11X2(3) ■ Alternatively the same signal could be searched 
for through the direct production of ncutralinos in asso- 
ciation with charginos in proton-proton interactions. 

This chain would give a clean final state signature [5(| 
with two opposite-sign same-flavour leptons. The kine- 
matics of the three-body decay imply that the invariant 
mass of any two particles in the decay must be smaller 
than the difference between the mass of the mother par- 
ticle and the mass of the remaining particle. Therefore, 
from the observation of one (or more) kinematic end- 
points in the invariant mass distribution of the two lep- 
tons from the decay given in Eq. (jTTJJ) the mass differ- 
ence m x q — m x <> can be measured. The error on this 
measurement has a systematic component, arising from 
lepton selection cuts, and a statistical component which 
will scale both with the value of the mass difference and 
with the number of events observed. This latter com- 
ponent will vary over the SUSY parameter space due to 
variations in the mass spectrum and in the total SUSY 
production cross-section. For the model studied in [49(, 
a precision of order 1% is obtained, which reflects the ex- 
cellent measurement capabilities of the LHC experiments 
for leptons. 

We have investigated the LHC reach for this measure- 
ment for the region of mSUGRA space defined in section 
4, by calculating the number of events expected at the 
LHC for the decay given in Eq. (fT0|) for an integrated 
luminosity of 100 fb _1 . The number of events shown are 
independent of the detector assumptions. For the quoted 
selection efficiencies we refer to an ATLAS analysis in 
parametrised simulation. The results are valid for the 
generic performance of an LHC detector and are expected 
to be equivalent for CMS. We used ISAJET 7.71 for 
the calculation of spectra and branching fractions and 
Prospino [57| for the NLO cross-section for gluino pro- 
duction. The results are shown in the right-hand side 
panel of Fig. [3] in the (mo — m i/2) plane before selection 
efficiency cuts. The sharp cut on the right-hand side is 
due to the fact that we do not consider models for which 
m x o — m x o > 85 GeV, at which point the kinematic end 
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m 1/2 (GeV) m 1/2 (GeV) 

FIG. 3: Left: Integrated 7-ray flux above 20 GeV, "I> 7 , in the (mo,nii/2) plane for the FP region defined in TableUwith Ao 
= 0, tan/3 = 10, sign(/i) = +1 and nit = 175 GeV. Overlaid is the H.E.S.S. II sensitivity for a 5cr detection in 50 hours 
(black contours). Right: Expected number of signal events before selection efficiency cuts for the leptonic 3-body decay of 
neutralinos (see Eq. (|10[) ) at the LHC for an integrated luminosity of 100 fb _1 in the (mo — mi/j) plane. 




FIG. 4: Left:Integrated 7-ray flux above 20 GeV in the FP region defined in Table U in the plane of the gluino mass, m§, as 
a function of the mass difference of the two lightest neutralinos m x o — m x o. Overlaid are the H.E.S.S. -II sensitivity (solid 
black contour) and the 5cr allowed range from WMAP QcdmIi 2 measurements (dashed purple contour). Right: Expected 
number of signal events before selection efficiency cuts for the leptonic 3-body decay of neutralinos (see Eq. HQ}) at the 
LHC for 100 fb _1 in the {m x o — m x o,nig) plane. 
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point in the two-lepton invariant mass distribution starts 
to be masked by the Z peak that arises from the decay 
in two leptons. When the decay \2 ^ ^Xi becomes 
kinematically accessible (i.e. when the neutralino mass 
difference exceeds the Z mass), the previous kinematic 
signature resulting from a three-body decay is lost, and 
the two body process that occurs instead does not al- 
low us to make an unambiguous and model independent 
measurement of the neutralino mass difference. This two 
body process will dominate up to the point where a neu- 
tralino decay channel featuring a light higgs h, \2 ~> ^Xii 
becomes open. The higgs decays mostly into bb pairs, 
and the extraction of the h — ► bb peak above the large 
combinatorial of b jets from gluino decays is very chal- 
lenging. For the three-body decay of the neutralinos to 
leptons, in the analysis of [49j an efficiency of ~ 25% is 
achieved for the signal, the SUSY background is ~ 0.6 
times the signal, and a SM background of ~ 250 events 
is expected for 100 fb _1 . The SM background is dom- 
inated by it production. This background, as well as 
the SUSY backgrounds, can be easily evaluated from the 
data themselves. In fact the two leptons from the three- 
body decays of the neutralino, being from the decay of 
a virtual Z, must have the same flavour and opposite 
sign. The two isolated leptons from it are produced from 
the decays of two W and can therefore have any flavour 
combination. It is therefore easy to statistically subtract 
the background by subtracting the invariant mass distri- 
bution of e^/Lt^ pairs from the sum of the distributions 
for e ± e =F and (i^n^ 1 pairs. Equivalent considerations are 
valid for the SUSY backgrounds, where the flavours of 
the two leptons are uncorrelated. Assuming that the ef- 
ficiency of the analysis cuts and ratio between signal and 
SUSY background do not vary dramatically in the inves- 
tigated region, a signal of ~350 events would be needed 
for the observation at 5cr level for an integrated luminos- 
ity of 100 fb _1 . The number of signal events is essen- 
tially determined by the gluino mass, as is clear in the 
right panel of Fig. [U where the number of signal events 
is given on the {m x o — m x o,jiij) plane. The opposite 
sign - same flavour two-lepton signal should therefore be 
observable for nig < 1300 GeV. 



VI. DISCUSSION OF THE 
COMPLEMENTARITY 

We discuss here the complementarity of the indirect 
detection with the phase II of H.E.S.S. and the ATLAS 
searches at LHC in the focus point region. 

The potentialities of the two experiments in the FP 
region can be readily compared through the integrated 
7-ray flux and the expected number of signal events for 
Next-to-LSP 3-body leptonic decays, in the (too,TOi/2) 
plane, displayed respectively in the left and right pan- 
els of Fig. [3J The left-hand side plot of Fig. [3] shows 
the H.E.S.S. 7-ray flux contours for a 5<r detection. The 
H.E.S.S. 7-ray integrated flux sensitivity which is at the 



level of ~ 10~ 12 cm _2 s _1 , allows us to test SUSY mod- 
els characterized by universal scalar soft mass parame- 
ters above 3300 GeV and gaugino soft mass parameters 
above 280 GeV. From the number of signal events for 
100 fb _1 in the 3-body leptonic decay shown in the right 
panel of Fig. [31 the ATLAS sensitivity will allow us to 
investigate mo lower than 4200 GeV and m 1 / 2 as high 
as 550 GeV. Thus the combination of the sensitivity of 
each experiment allows us to cover the overall region. 
This point highlights the complementarity between the 
indirect detection of dark matter and the SUSY searches 
at LHC. Moreover, an overlap is obtained in a small re- 
gion characterized by scalar soft masses in the 3300-4200 
GeV range and gaugino soft masses between 280 and 550 
GeV. In case of a signal discovery consistent with this re- 
gion, complementary cross-checks between the two exper- 
iments would lead to further constraints on the allowed 
parameter space in the (too,TOx/2) plane. 

For further illustration, we have chosen in Fig. [H a 
given point in the parameter space: 

m = 3350 GeV, m 1/2 = 300 GeV, A = GeV, 

/i > 0, tan/3 = 10 (11) 

This point, with the top mass set to 175 GeV and the 
mass spectrum computed with ISAJET, is close to the 
point (sometimes dubbed SU(2) ) which has been chosen 
by the ATLAS experiment for a detailed study [49j . For 
each quantity reported in Table 11111 two values are given 
corresponding to the output spectra of ISAJET7.69 and 
ISAJET7.71. This allows to give an estimate for the 
theoretical uncertainties. As can be seen from Fig. [Hand 
Table [Hll this particular point, which is within the reach 
of ATLAS, is also within the reach of H.E.S.S. II provided 
a typical 50 hours observation time of a Galactic Center 
type source. 

As explained in section [Vl a key parameter in the AT- 
LAS sensitivity is the gluino mass. The sensitivity curves 
of each experiment allow to constrain different regions in 
the plane of the gluino mass versus the mass difference 
of the two lightest neutralinos as shown in Fig. [U As can 
be seen, one anticipates lower flux values for larger neu- 
tralino mass differences. Given an astrophysical model 
for the source (see section 4), H.E.S.S. II would thus 
put lower limits on the neutralino mass differences which 
can be further refined with some knowledge of the gluino 
mass. Moreover, given the sensitivity zone of H.E.S.S. 
II, one finds that a larger gluino mass would allow the 
exploration of a larger fraction of the parameter space. 
In contrast, and as illustrated in from Fig.HJ ATLAS will 
have larger sensitivities for lighter gluinos. In any case, 
and irrespective of the gluino mass, the extraction of in- 
formation from ATLAS measurements becomes problem- 
atic for a neutralino mass difference exceeding 85 GeV. 
Finally, comparing the two considered observables in the 
panels of Fig. [4] one notes a reversed sensitivity to the 
mass parameters: ATLAS would provide a strong con- 
straint on the gluino mass whereas the indirect detection 
is more sensitive to the neutralino mass differences. Or- 
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TABLE III: Values of the neutralino relic density Q x h 2 , the 7-ray flux $ 7 , the number of events per 100 fb -1 after selection 
efficiency cuts for NLSP 3-body leptonic decays, the mass difference of the two lightest neutralinos and the gluino mass, for the 
SUSY point defined in Eq. (jTTJ) . The two quoted values for each observable correspond respectively to the output of IS A JET 
7.69 and 7.71 (see text for details). 



Relic density Q. x h 2 


7-ray flux $ 7 (cm 2 s x ) 


Events per 100 fb" 1 


m^o — m<a (GeV) 


m- g (GeV) 


0.06 0.16 


6.6xl0~ 12 2.3xl0 -12 


855 468 


55 88 


852 



thogonal constraints can thus be obtained by indirect 
detection and collider searches, and their possible combi- 
nation would reduce significantly the allowed parameter 
space. 

VII. CONCLUSION AND PERSPECTIVES 

In this paper we have examined the complementar- 
ity of two approaches to unravel signatures of dark mat- 
ter particles in the case where the latter are neutralino 
LSPs in the Focus Point region of the mSUGRA sce- 
nario. Searches at colliders in this region which is dif- 
ficult to investigate as compared to the so-called bulk 
region, will benefit from complementary input from indi- 
rect searches through 7-rays|60j. In the present study we 
used the most realistic simulations for both ATLAS and 
H.E.S.S. II for which we considered the Galactic Center 
as a benchmark astrophysical DM source. We found that 
the two experiments will typically explore two different 
regions of the studied parameter space with some possi- 
ble overlap, allowing in principle complementary analy- 
ses. We illustrated these features in terms of the universal 
soft parameters as well as in terms of the physical masses, 
the first corresponding to a top-down model-dependent 
approach while the second illustrates the more challeng- 
ing model-independent one. For instance, in the scanned 
range of neutralino LSP masses from 50 to 300 GeV, AT- 
LAS will be sensitive to gluinos lighter than 1300 GeV 
whereas H.E.S.S. II, with a flux sensitivity of the order of 



10~ 12 cm _2 s _1 , will be able to cover all the region above 
1 TeV. Furthermore, the mass difference of the two light- 
est neutralinos which is a key observable for ATLAS, is 
interestingly found to be a sensitive parameter for the 
7-ray fluxes that H.E.S.S. II can observe or constrain. 

More generally, in the very near future, H.E.S.S. II and 
LHC will be two major experiments allowing us to probe 
the supersymmetric hypothesis for dark matter with an, 
up to now, unequalled capability. Neutralino masses in 
the yet uncovered 100 GeV range will be accessible to 
H.E.S.S. II thus filling the energy gap between the cur- 
rent H.E.S.S. I and the GLAST satellite experiment to 
be launched in Spring 2008. Finally, if SUSY signals are 
discovered at the LHC, the cross-sections, branching ra- 
tios and masses of the new particles determined with a 
few fb -1 luminosity, would provide valuable search win- 
dows both for direct and indirect supersymmetric dark 
matter detection. Of course, substantial systematic un- 
certainties arise on the astrophysical parameters in indi- 
rect searches, which can be reduced through an improved 
knowledge of the halo profiles of the observed sources. 
Direct detection also suffers from systematic effects re- 
lated to astrophysical and nuclear parameter uncertain- 
ties which affect the constraints on the nucleon-WIMP 
cross-sections. The more model-independent measure- 
ments of the SUSY parameters will require data samples 
collected with hundreds of fb -1 luminosity at the LHC, 
as well as independent input from indirect (or direct) 
dark matter searches as illustrated in this paper. 
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the acronym mSUGRA is used everywhere in this paper 
in a loose sense, that is without assuming a GUT scale re- 
lation between the soft parameters Aq and Bo that would 
result from a flat Kahler metric. 

More precisely the most important differences occur 
in fact between Isajet 7.69 on one side and Sus- 
pect/SoftSuSy Q/Spheno on the other for the de- 
termination of the fj, parameter, while the differences 
among the latter three codes appear to be generally much 
milder [41J . It should be noted however that the discrep- 
ancies in the LSP neutralino mass remain moderate even 
though in the focus point region its higgsino component 
is not negligible, since the neutralino N\ component re- 
mains essentially determined by the U(1)y gaugino soft 
mass term Mi in this region, so that its mass is not very 
sensitive to fj,. 

We note that searches in the bulk region are easier due to 
the fact that the lower value of mo ensures that squarks 
and sleptons will be visible at the LHC. The greater range 
of visible sparticles will provide a greater knowledge of 
the sparticle mass spectrum. 



